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The  inclusion  rating  method  by  statis- 
tics  of  extreme  values  (IRMSE)  using 
y/area  of  inclusions  as  the  size  parame- 
ter enables  one  to  discriminate  between 
current  super-clean  steels.  Moreover, 
IRMSE  enables  one  to  predict  the  size 
(\/areami,x)  of  maximum  inclusions  con- 
tained in  domains  larger  than  the  in- 
spection domain.  The  statistical 
distribution  of  y/aream^x  can  be  used 
for  the  quality  control  of  materials  and 


for  the  prediction  of  a scatter  band  of 
fatigue  strength.  Practical  procedures  of 
inclusion  rating  and  prediction  of  a 
scatter  band  of  fatigue  strength  are 
shown. 
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1.  Introduction 

With  the  increase  in  cleanliness  of  steels,  con- 
ventional inclusion  rating  methods  are  no  longer  as 
useful  as  before,  because  conventional  inclusion 
rating  methods  cannot  determine  the  cleanliness  of 
new  clean  steels.  Although  the  cleanliness  of  steels 
has  been  markedly  improved  in  the  last  two 
decades,  the  fatigue  strength  of  recent  clean  high 
strength  steels  cannot  attain  the  ideal  value  ex- 
pected from  their  high  static  strength.  Nonmetallic 
inclusions  are  predominantly  the  cause  of  lower  fa- 
tigue strength  even  for  such  clean  high  strength 
steels.  Thus,  in  order  to  predict  the  fatigue  strength 
behavior  and  to  evaluate  quality,  we  need  a new 
inclusion  rating  method  relevant  to  recent  super- 
clean steels.  The  inclusion  rating  method  based  on 
statistics  of  extreme  [1]  is  most  relevant  for  this 
purpose.  In  the  following,  we  call  this  method  In- 
clusion Rating  Method  by  Statistics  of  Extreme 
(IRMSE). 

In  this  study,  we  shall  first  show  that  if  we  choose 
an  appropriate  size  parameter  for  inclusions,  the 


size  of  inclusions  obey  the  statistics  of  extreme 
value  theory.  The  appropriate  size  parameter  is  the 
square  root  of  projected  area  of  the  maximum  in- 
clusion contained  in  a standard  inspection  area  or 
volume,  \/aream»x.  Second,  we  predict  the  size  of 
the  maximum  inclusion  which  may  be  contained  in 
a larger  area  or  volume  than  the  standard  inspec- 
tion  area  and,  lastly,  we  use  the  size  parameter, 
\/areamax,  to  predict  the  scatter  band  of  fatigue 
strength  of  hard  steels. 

The  merits  of  IRMSE,  in  comparison  with  con- 
ventional methods,  are  (1)  to  distinctly  discrimi- 
nate the  cleanliness  of  recent  super-clean  steels, 
and  (2)  to  predict  the  size  of  larger  inclusions  con- 
tained in  a domain  larger  than  the  inspection  do- 
main. This  method  is  useful  for  quality  control  of 
materials  and  for  improvement  of  the  steel  making 
processes.  It  also  enables  one  to  predict  the  scatter 
of  the  fatigue  strength  of  a large  number  of  mass 
production  products. 
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2.  Nonmetallic  Inclusions  as  a Fatigue 
Fracture  Origin 

Figure  1 shows  an  example  of  the  nonmetallic 
inclusion  which  was  observed  at  fatigue  origin  of  a 
bearing  steel  under  a rotating  bending  fatigue  test. 
If  this  inclusion  did  not  exist  in  this  specimen,  the 
fatigue  strength  of  this  specimen  should  have  been 
higher  than  the  applied  stress,  cr„  = 1078  MPa. 
Since  the  size  and  location  of  nonmetallic  inclu- 
sions scatter  randomly,  the  fatigue  strength  of  high 
strength  steels  naturally  scatters.  Although  there 
has  been  a firm  opinion  that  the  chemical  composi- 
tion and  shape  of  nonmetallic  inclusions  substan- 
tially influences  the  fatigue  limit,  Murakami  et  al. 
[2-5]  have  shown  the  incorrectness  of  the  conven- 
tional opinion  by  their  detailed  experiments  and 
analyses,  and  reported  distinct  experimental  evi- 
dent that  the  size  of  inclusions  (defined  by 
\/area ) is  the  most  crucial  geometrical  parameter. 
It  is  empirically  known  that  the  intrinsic  fatigue 
strength  of  steels  is  determined  by  the  hardness 
(Hv)  of  its  microstructure.  For  steels  with  //v<400, 
nonmetallic  inclusions  contained  in  current  com- 
mercial steels  are  not  detrimental  and  we  have  the 
following  empirical  formula 

(t^  = 1.6Hv  (1) 


Vickers  hardness  //v  = 745  kgf/mm’ 

Applied  stress  at  surface  r/o  = 1078  MPa 
Cycles  to  failure  Nf  = 7.94x  10'’ 

Square  root  of  projection  area  of  inclusion  \/area  =9.8  jcm 
Distance  from  surface  /i  = 158  jim 
Applied  stress  at  inclusion  a-'  = 1034  MPa 
Chemical  composition  of  inclusion:  Al-Mn-S-O 


where  a^,  is  the  fatigue  limit  (MPa)  and  Hv  is  the 
Vickers  hardness  (kgf/mm^).  However,  for  steels 
with  //v>400,  the  effect  of  inclusions  reveals  itself 
and  the  intrinsic  or  ideal  fatigue  limit  given  by  Eq. 
(1)  cannot  be  attained.  The  fatigue  strength  de- 
pends on  the  size  {\/area ) and  location  of  the  fatal 
inclusion  and  Hv  of  the  matrix.  Murakami  et  al.’s 
[6-9]  fatigue  limit  prediction  equations  are  classi- 
fied into  three  categories  depending  on  the  loca- 
tion of  fatal  inclusions  (see  Fig.  2): 

Fatigue  limit  for  a surface  inclusion  [Fig.  2(a)] 

tr.  = 1.43  {Hv+m)l{V^y'^  (2) 

Fatigue  limit  for  an  inclusion  in  touch  with  free 
surface  [Fig.  2(b)] 


a-„,  = 1.41  {Hv  + myiVarea  )‘'^  (3) 

Fatigue  limit  for  an  internal  inclusion  [Fig.  2(c)] 


O'.  = 1.56  (Hv  + m)f{Varea  )‘^^  (4) 

where  the  units  are  cr»,:  MPa,  \/area:  |xm,  and  Hv'- 
kgf/mm^. 


Free  surface  Free  surface  Free  surface 


(a)  Surface  inclusion  (b)  Inclusion  just  (c)  Interior 

below  surface  inclusion 

Fig.  2.  Various  locations  of  inclusions  causing  fatigue  fracture. 

Since  for  a constant  value  of  area,  an  inclusion  is 
most  detrimental  when  it  exists  just  in  touch  with 
the  free  surface  of  a specimen,  we  can  use  Eq.  (3) 
in  combination  with  the  maximum  size  VflfCflmax 
obtained  by  IRMSE  to  predict  the  lower  bound 
(cr„,^)  of  scattered  fatigue  strength  of  many  speci- 
mens or  machine  elements. 


3.  Inclusion  Rating  of  Various  High 
Strength  Steels  by  Statistics  of 
Extreme 

Figure  3 explains  the  practical  procedure  to  im- 
plement the  inclusion  rating  by  statistics  of  extreme 
values.  The  details  of  this  method  are  reported  in 


Fig.  1.  A typical  example  of  inclusion  observed  at  the  center  of 
fatigue  fracture  origin  [super-clean  bearing  steel,  SUJ2(N)]. 
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Fig.  3.  Practical  procedure  of  the  inclusion  rating  by  statistics  of 
extreme  values. 

Murakami  et  al.’s  papers  [3,5,9-11].  The  procedure 
is  briefly  explained  in  the  following,  see  [11]. 

(1)  A section  perpendicular  to  the  maximum 
principal  stress  is  cut  from  the  specimen.  After  pol- 
ishing with  a n°2000  emery  paper,  the  test  surface 
is  mirror-finished  with  buff. 

(2)  A standard  inspection  area  5o  (mm^)  is  fixed. 

Generally,  it  is  advisable  to  take  a microscope  pic- 
ture for  an  area  approximately  equivalent  to  So.  In 
the  area  So,  the  inclusion  of  maximum  size  is  se- 
lected Then,  the  square  root  of  the  projected  area 
y/areomix  of  this  selected  inclusion  is  calculated. 
This  operation  is  repeated  n times  (in  n areas  So) 
(see  Fig.  2).  

(3)  The  values  of  y/aream^j  are  classified,  start- 
ing from  the  smallest,  and  indexed:  (with  j = 

We  then  have  the  following  relation: 


y/area^„,i^  \/«^^«max.2^  y/area:„^,„. 

The  cumulative  distribution  function  Fj  and  the 
reduced  variates  y,  are  then  calculated  from  the 
equations. 

Fj  =j  X 100/(n  + 1) 
yj  = - ln[  - ln(/7(/i  -bl)]. 

(4)  The  data  are  then  plotted  on  probability  pa- 
per.  The  point  j has  an  abscissa  coordinate  of 
y/areumaxj  while  the  ordinate  axis  represents  either 
Fj  oryj.  An  example  of  the  curve  is  shown  in  Fig.  4. 

Figure  4 shows  the  inclusion  ratings  by  IRMSE 
for  two  kinds  of  super-clean  bearing  steels, 
SUJ2(N)  and  SUJ2  (H).  The  total  oxygen  con- 
tained in  these  steels  is  8 ppm  for  SUJ2(N)  and 
5 ppm  for  SUJ2(H).  This  kind  of  information 


Fig.  4.  Cumulative  frequency  of  the  extreme  values  of  inclusions 
[Super-clean  bearing  steels,  SUJ2(N)  and  SUJ2(H)]. 

enables  one  to  discriminate  quantitatively  the  dif- 
ference among  the  cleanliness  levels  of  the  same 
kind  of  materials  produced  by  different  companies 
or  produced  by  a company  at  different  periods. 
Thus,  this  information  will  be  useful  for  the  quality 
control  of  materials  and  the  improvement  of  the 
steel  making  process. 

It  is  not  a priori  evident  to  what  extent  the  ex- 
treme values  y/areum^  of  inclusions  contained  in 
various  steels  follow  extreme  statistics  value.  How- 
ever, Murakami  et  al.  [3,5-11]  have  shown  many 
examples  of  measurements  which  obey  the  statis- 
tics of  extreme  value  theory.  Uemura  and  Mu- 
rakami [12]  carried  out  a three-dimensional 
numerical  simulation  to  find  the  statistical  distribu- 
tion of  the  extreme  values  y/aream^  of  inclusions 
which  were  distributed  in  a constant  volume  with 
the  size  (D)  distribution  of  the  type, 

0(D)  =—  exp  (-— ), 
m m 

where  m is  the  mean  value,  and  they  confirmed  the 
validity  of  IRMSE  (Fig.  5).  In  addition,  they  indi- 
cated the  quantitative  difference  between  two-di- 
mensional and  three-dimensional  measurements, 
though  the  difference  virtually  vanishes  with  in- 
creasing inspection  domains. 


347 


Volume  99,  Number  4,  July-August  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


. . ■ ■ 1 200/mi  I 

(a-1)  Inclusion  distribution  on 
inspection  section  of  the 
material  with  m=l  ^m. 


. 1 20b>mi  |- 

(b-1)  Inclusion  distribution  on 
inspection  section  of  the 
material  with  m=2  ^.m. 


. • .•  1 200//rn  | 

(c-1)  Inclusion  distribution  on 
inspection  section  of  the 
material  with  m=3  ^m. 


F(%) 


(a-2)  Statistical  distribution  of 
true(3-D)  and  apparent 2-D) 
extreme  values  of  inclusions 
for  material  with  m=l  [^m 
(So=0.25  mm2,  number  of 
inspections  = 40). 

F(%) 


(b-2)  Statistical  distribution  of 

true(3-D)  and  apparent(2-D) 
extreme  values  of  inclusions 
for  material  with  m=2  (xm 
(So=0.25  mm2,  number  of 
inspections  = 40). 


F(%) 


(c-2)  Statistical  distribution  of 

true(3-D)  and  apparent(2-D) 
extreme  values  of  inclusions 
for  material  with  m=3  jxm 
(S(j=0.25  mm2,  number  of 
inspections  = 40). 


Fig.  5.  Numerical  simulation  of  the  inclusion  rating  by  statistics  of  extreme  values  on  the  materi- 
. . 1 D 

als  with  the  inclusion  size  distribution  of  the  type  4>(D)= — exp  ( ). 

m m 
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4.  Application  to  Prediction  of  Scatter 
Band  of  Fatigue  Strength 

Figure  6 illustrates  the  shape  and  dimension  of  a 
tension-compression  fatigue  specimen  [13].  The 
material  used  is  tool  steel,  SKH51.  The  chemical 
composition  is  shown  in  Table  1.  Table  2 shows  the 
mechanical  properties. 


Fig.  7.  Statistical  distribution  of  the  extreme  values,  the  maxi- 
mum size  of  inclusion  at  the  center  of  fracture  origin  (Tool 
steel,  SKH51). 


Fig.  6.  Shape  and  dimension  of  tension-compression  fatigue 
specimen  (mm)  (Tool  steel,  SKH51). 

Figure  7 shows  the  extreme  value  distribution  of 
y/area  of  the  inclusions  found  at  the  fracture  origin 
of  34  specimens.  The  data  in  Fig.  7 are  the  extreme 
values  obtained  by  the  fatigue  test  but  not  by  the 
two-dimensional  metallographic  method  described 
in  Sec.  3.  Figure  8 indicates  the  location  of  these 
inclusions  on  the  fracture  surface.  If  the  tension- 
compression  fatigue  test  is  not  performed  correctly, 
that  is,  specimens  are  subject  to  a bending  moment 
due  to  a bad  alignment  or  the  curving  of  the  speci- 
men axis,  nonmetallic  inclusions  existing  near  the 
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Table  1.  Chemical  composition  in  wt%  of  material  (Tool  Steel, 
SKH51) 


c 

Si 

Mn 

P 

s 

Cr 

W 

Mo 

0.81 

0.31 

0.29 

0.018 

0.002 

3.92 

6.10 

4.85 

V 

Co 

Cu 

Ca 

A1 

Mg 

O 

1.81 

0.46 

0.07 

0.004 

0.035 

0.0005 

0.0018 

Fig.  8.  Relationship  between  the  size  (y/area)  and  loeation  of 
inclusions  at  the  center  of  fraeture  origin  (Tool  steel,  SKH51). 

free  surface  are  likely  to  appear  as  the  fracture 
origin  on  the  fracture  surface  [14].  In  such  a case, 
unusually  low  fatigue  strength  is  likely  to  be  ob- 
tained. Since  the  fatigue  fracture  origins  shown  in 
Fig.  8 are  distributed  randomly  on  the  section  of 
specimen,  these  data  may  be  valid  for  the  statistical 
analysis.  However,  it  should  be  noted  that  when 


Table  2.  Mechanical  properties  of  quenched  and  tempered  test  material  (Tool  steel, 
SKH51) 


Vickers 

Heat 

0.2% 

Tensile 

Reduction 

hardness 

treatment 

Proof  stress 

strength 

Elongation 

of  area 

//v 

(MPa) 

(MPa) 

(%) 

(%) 

(kgf/mm^) 

Heat  treat.  1 

1820 

2 no 

2.9 

3.7 

615 

Heat  treat.  3 

2 270 

2 560 

2.0 

0 

654 
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the  surface  inclusions  became  the  fracture  origins, 
the  data  were  not  plotted  on  Fig.  7,  because  such 
inclusions  are  more  detrimental  than  an  inclusion 
having  the  same  size  and  existing  internally  and  ac- 
cordingly they  may  be  a little  smaller  than  the  exact 
maximum  inclusion. 

In  the  case  of  the  data  of  Fig.  7,  the  volume  of 
the  test  part  of  one  specimen  (Fig.  6)  corresponds 
to  one  inspection  domain  and  there  are  34  extreme 
values  in  Fig.  7.  Therefore,  Fig.  7 can  be  used  for 
predicting  the  expected  maximum  size  of  the  inclu- 
sion which  may  be  contained  in  more  specimens 
than  those  used  in  fatigue  tests.  For  example,  an 
inclusion  having  \/^''^^^max=138  jo-m  is  expected  to 

be  contained  in 100  specimens  (N  = 100). 

Combining  this  \'aream^  ( = 138  ixm)  and  Eq.  (3), 
the  lower  bound  {cr^e)  of  fatigue  strength  of  100 
specimens  can  be  predicted. 

Figure  9 compares  the  scatter  observed  in  exper- 
iments and  the  predicted  lower  bound  (r„e  of  the 
scatter  band.  The  prediction  is  in  good  agreement 
with  experiments.  The  prediction  of  the  lower 
bound  of  fatigue  strength  explained  above  can  be 
used  for  the  quality  control  of  machine  elements 
which  are  produced  by  mass-production  and  can- 
not be  tested  individually. 

The  data  as  shown  in  Fig.  7 offer  us  reliable  in- 
formation on  inclusions  expected  to  be  contained 
in  other  specimens.  However,  obtaining  the  data 
shown  in  Fig.  7 requires  preparation  of  many  pre- 
cise specimens  and  time  consuming  fatigue  tests. 
To  avoid  this  inconvenience,  the  author  has  pro- 
posed an  alternative  two-dimensional  method  as 
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Fig.  9.  Comparison  between  the  experimental  results  and  the 
lower  bound  of  fatigue  strength  which  was  predicted  on  the  ba- 
sis of  Eq.  (3)  and  the  maximum  size  of  inclusion  (Tool  steel, 
SKH51). 


explained  in  Sec.  3.  A sufficient  number  (TV)  of  in- 
spection domains  (inspection  areas)  necessary  to 
predict  reliably  for  more  specimens  or 

larger  areas  should  depend  on  the  materials  to  be 
inspected  and  on  the  inspection  area  So  observed 
by  the  image  processor  combined  with  an  optical 
microscope.  From  the  author’s  experience,  it  is  rec- 
ommended that  TV  be  larger  than  40  for  5o  = 0.031 
mm^. 

Several  Japanese  industries  have  already  put  the 
method  proposed  in  this  study  in  practice  [15]. 


5.  Conclusions 

(1)  If  we  define  the  size  of  nonmetallic  inclu- 
sions contained  in  commercial  steels  by  the  square 
root  of  the  projected  area,  \/area,  the  maximum 
values,  \/aream^,  in  a definite  inspection  domain 
obey  the  statistics  of  extreme  value  theory. 

(2)  The  inclusion  rating  method  by  the  statistics 

of  extreme  values  (IRMSE)  based  on  can 

be  used  for  a new  inclusion  rating  method.  IRMSE 
enables  one  to  discriminate  distinctly  between  re- 
cent super-clean  steels,  while  conventional  inclu- 
sion rating  methods  are  no  longer  valid  as  the 
method  to  evaluate  the  cleanliness  of  new  clean 
steels. 

(3)  IRMSE  is  useful  not  only  for  a relative  eval- 
uation of  materials  but  also  for  the  prediction  of 
the  expected  maximum  size  of  inclusions  to  be  con- 
tained in  a domain  larger  than  the  inspection  do- 
main. The  value  of  \/areamM  can  be  used  with  the 
fatigue  strength  prediction  equation  to  predict  a 
scatter  band  of  fatigue  strength  of  high  strength 
steels. 
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